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We consider an evolution operator for a discrete Langevin 
equation with a strongly hyperbolic classical dynamics and 
noise with finite moments. Using a perturbative expansion 
of the evolution operator we calculate high order corrections 
to its trace in the case of a quartic map and Gaussian noise. 
The leading contributions come from the period one orbits 
of the map. The asymptotic behaviour is investigated and is 
found to be independent up to a multiplicative constant of 
the distribution of noise. 



I. INTRODUCTION 

In natural phenomena stochastic processes of various 
strength always have an influence. In the three preceding 
papers the effects of noise on measurable properties 
such as dynamical averages Q in classical chaotic dy- 
namical systems were systematically accounted. In order 
to construct a statistical theory of dynamical systems 
one introduces densities of particles governed by a corre- 
sponding evolution operator. For a repeller the leading 
eigenvalue of this operator C yields a physically measur- 
able property of the dynamical system, the escape rate 
from the repeller. In the case of deterministic flows, the 
periodic orbit theory yields explicit and numerically effi- 
cient formulas for the spectrum of C as zeros of its spec- 
tral determinant 

The theory developed is closely related to the semi- 
classical expansions based on Gutzwiller's formula for 
the trace in terms of classical periodic orbits ||] in that 
both are perturbative theories (in the noise strength or 
Ti) derived from saddle point expansions of a path inte- 
gral containing a dense set of unstable stationary points 
(typically periodic orbits). The analogy with quantum 
mechanics and field theory is made explicit in |y| where 
Feynman diagrams are used to find the lowest nontrivial 
noise corrections to the escape rate. 

An elegant method, inspired by the classical pertur- 
bation theory of celestial mechanics, is that of smooth 
conjugations In this approach the neighbourhood of 
each saddle point is flattened by an appropriate coordi- 



nate transformation, so the focus shifts from the orig- 
inal dynamics to the properties of the transformations 
involved. The expressions obtained for perturbative cor- 
rections in this approach are much simpler than those 
found from the equivalent Feynman diagrams. Using 
these techniques, we were able to extend the stochas- 
tic perturbation theory to the fourth order in the noise 
strength. 

In 1^] we develop a third approach, based on construc- 
tion of an explicit matrix representation of the stochas- 
tic evolution operator. The numerical implementation 
requires a truncation to finite dimensional matrices, and 
is less elegant than the smooth conjugation method, but 
makes it possible to reach up to order eight in expansion 
orders. As with the previous formulations, it retains the 
periodic orbit structure, thus inheriting valuable infor- 
mation about the dynamics. 

The corrections to the escape rate were found to be a 
divergent series in the noise expansion parameter. This 
reflects that the corrections are calculated using the so- 
called cumulant expansion from other divergent quanti- 
ties, the traces of the evolution operator £" j^]. In this 
article the focus is on the high order corrections for the 
special case of the first trace, tr(£). For the asymptotic 
study we have to calculate a sufficient number of correc- 
tions and here we are led naturally to a contour integral 
method. Next the asymptotic behaviour is extracted by 
the method of steepest descent. 

II. THE STOCHASTIC EVOLUTION OPERATOR 

In this section we introduce the noisy repeller and its 
evolution operator. 

An individual trajectory in presence of additive noise 
is generated by iterating 

Xn+l = f{Xn) + Crin , (1) 

where f{x) is a map, ^„ a random variable with the nor- 
malized distribution p(<^), and a parametrizes the noise 
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strength. In what follows we shall assume that the map- 
ping f{x) is one-dimensional and expanding, and that 
the ^„ are uncorrelated. A density of trajectories 0(a;) 
evolves with time on the average as 

0n+l(y) = (-Co^n) (y) = J dxC{y,x)(t)n{x) (2) 

where the £ evolution operator has the general form of 



C{y,x) = 6„{y - f{x)), 



(3) 



For the calculations in this paper Gaussian weak noise is 
assumed. In the perturbative limit, cr 0, the evolution 
operator becomes 
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where (5(2") denotes the 2n-th derivative of the delta dis- 
tribution and a2n — (2"n!)~^ or a2„ — -^f^ for general 
noise with finite n'th moment to„. Strictly speaking we 
here only consider distributions with vanishing odd mo- 
ments, but the formulas below can easily be generalized. 
The perturbative form of the evolution operator is ob- 
tained formally from (j^) by Taylor expansion or by doing 
a saddle point integral discarding subdominant terms. In 
particular there may exist stationary points which do not 
correspond to periodic orbits. However, these will give 
exponentially small contributions and are therefore omit- 
ted in the weak noise limit cr ^ 0. The map considered 
here is the same as in our previous papers, a quartic map 
on the (0, 1) interval given by 



f{x) = 20 



1 
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(6) 



we reformulate expression as 



= lim y 



(2n)!a„(T2 



2m 



dx 
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dx 



[x - f{x) + ie)2"+ 



[x - f{x) - ie)2»+i 

(9) 

section 5.3 equation 



obtaining just the result as in 
(19) p.119. 

The integrals in (^) have poles close to the fixpoints of 
the map, which are defined by 



fix*)^0. 



(10) 



For the map in question, a quartic, there are four points 
in the complex plane which satisfy this equation: two 
of them are on the real axis: one at xq = 0, the other 
xi close to unity, as it can be seen on figure (0). The 
two real fix points correspond to orbits and 1 in the 
language of symbolic dynamics . 




FIG. 1. The map on the [0,1] interval 



The integral in (|^) containing ~ie in the denominator 
has a pole at a;o shifted down from the real axis into the 
negative imaginary half plane and a pole at xi shifted 
up into the positive imaginary half plane. At the other 
integral in (^ the poles are moved in the opposite di- 
rection. We shall now evaluate the integrals by contour 
integration. The chosen contour is shown on figure (0). 



III. RESIDUE METHOD FOR THE TRACE 

The trace of the evolution operator in the weak-noise 
limit is obtained from 



/OO 
da;^a2„cr2"(5(2")(a;_/(a;)). 

n=0 

By expressing the delta function as 
6{x) — lim — Im 



(7) 
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(8) 
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FIG. 2. The chosen contour and the motion of poles. 

The integral containing the —ie will pick up contribu- 
tion from the fixpoint at x — 0, the other integral from 
the fixpoint at xi. Both will pick up a contribution from 
the complex fixpoint, but they cancel out in the e — > 
limit. The result of the integration is 



where 



n=0 



(2n)!cr' 
2«n! 



[Resi^„(a:o)-Resi^„(a;i)], (11) 
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(12) 



We used this method to calculate further corrections 
(around cr™) to the trace for the orbit and orbit 1. This 
exact method justified the previous approximate calcu- 
lations to all known digits. It turns out that the cor- 
rections for orbit all were positive whereas those for 
orbit 1 involved sign changes with an apparent period. 
Furthermore both were asymptotic series. 



IV. HIGHER ORDER NOISE CORRECTIONS 

Practical reasons restricted the calculations with the 
residue method to around order seventy in the expansion 
parameter. Nevertheless it is possible to give approxi- 
mate statements for even higher order. 

First we transform the residue integrals in (||) as 



(2?T,)!a2n g 

27rr 



2n 



dx 



{x - fix) T ie)2"+i 



(2n)!a2„Cr^" I ^^^{~(2n+l)log(x-f{x)Tte)) 

2ni J 



(13) 



Using n as a large parameter the integral is next cal- 
culated by the method of steepest descent. The initial 
loops, Co and Ci, for each fixpoint are expanded until the 
paths of steepest descents are reached. We refer to fig. ||. 
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FIG. 3. Contourlines of \x — ,f{x)\. Filled circles indicate 
fixpoints of the map, crosses correspond to saddles and thick 
lines show paths of steepest descent. 

The paths of steepest descent pass through the saddle 
points satisfying 

/■(a,,) = 1. (14) 

Around each saddle we calculate the leading contri- 
bution and its corrections as usual. In the following we 
will divide the discussion into orbit and orbit 1. For 
our map there are three saddle points: one real and two 
complex conjugate. 

In case of orbit 0, the integral is carried out on the 
infinite path Vq shown on figure fig. ^. Here only the 
saddle point on the real axis contributes. 

The result of the integration to leading order is 



2"+in!7r 



{2n+l)f"{xs) 



2"r(n + 1/2) 



2nixs - /(x.))2"^2(n + l/2)(x, - /(a;,))/"(a;,) 

(15) 

For a convenient representation we introduce the factors 

1 



Co 



2TT^2iXs~fiXs))f'{Xs)y 

2 



\x-f{x)Y 
and the result takes the form 

cr2"6Jfr(n-|- 1/2) 



v/n+ 1/2 



(16) 
(17) 

(18) 



We show the ratio of consecutive terms /„_|_i//„ on 
fig. with the ratios from the exact residue calculation 
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and conclude that the real saddle indeed controls the 
leading behaviour. 



FIG. 4. Ratio of noise corrections I^?-tl for orbit 0. 

-'211 

For the orbit 1 all three saddles are relevant, since the 
initial loop Ci decomposes into Vq^Vi and 7^2- However, 
the real saddle now becomes subdominant as can be seen 
from e.g. the contour plot fig. |[ By analyticity of the 
map the two complex saddles give results that are com- 
plex conjugate to each other. 

At the complex saddle points (zi, Z2) we introduce sim- 
ilar factors Ci,6i as in (p^,p^) now complex. 

The contribution from these two complex saddles to 
can now be written as 



r(n+l/2) 



r{n- 



-1/2 
1/2) 



^/n + 1/2 



[Ci6? + C+(6+)"] 
[rq" cos(n0 — a)] . 



(19) 



Here r, g, and a are easily found from Ci, bi above. 

The final result for orbit 1 also includes the real saddle 
and to leading order is given by 



r(n + 1/2) 
v/n+1/2 



Co^o + cos{n(j) - a)] 



(20) 



We remark that the same calculation could be done for 
another distribution than Gaussian. One would find an 
identical functional form up to a multiplicative constant, 
the corresponding moment. 




FIG. 5. Ratio of noise corrections 



for orbit 1 



Here the plot of the ratios of the consecutive terms 
fig. 1^ show a more complex behaviour than for orbit 0. 
This we attribute to the presence of dominant complex 
saddles, which will introduce phases and hence the sign 
changes already mentioned. The apparent periodicity of 
seven is nicely captured by the two complex saddles. Nev- 
ertheless, high order calculations actually find deviations 
from this pattern in agreement with that the phase cal- 
culated does not evaluate exactly to an integer multiple 
of 2tt/7. The saddle point expansion in general improves 
for high orders as can be seen on fig. o. For our calcu- 
lations two corrections to the complex saddles give the 
best results in the regime of n we are considering. This 
is because as an asymptotic series truncation is required 
at an optimal point when n is finite. So adding more 
corrections will only improve the results for higher n but 
actually make the lower less accurate. 
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FIG. 6. The logarithm of relative error for the noise cor- 
rections for orbit 1 



4 



V. SUMMARY AND CONCLUSION 



We have studied the evolution operator for a discrete 
Langevin equation with a strongly hyperbolic classical 
dynamics and noise with finite moments. Using a per- 
turbative expansion of the evolution operator we have 
calculated high order corrections to its trace in the case 
of a quartic map and Gaussian noise. The leading contri- 
butions come from the period one orbits of the map. We 
have found the asymptotic behaviour of the corrections 
using the method of steepest descent. The functional 
form of the high order corrections have the same form 
independent of the actual noise distribution up to a con- 
stant given by the corresponding moment. The asymp- 
totics of the trace of the evolution operator are governed 
by subdominant terms corresponding to terms previously 
neglected in the perturbative expansion. 
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